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Free nucleons propagating in water are known to produce 7 rays, which form a background 
to the searches for diffuse supernova neutrinos and sterile neutrinos carried out with Cherenkov 
detectors. As a consequence, the process of nucleon knockout induced by neutral-current quasielastic 
interactions of atmospheric (anti)neutrinos with oxygen needs to be under control at the quantitative 
level in the background simulations of the ongoing and future experiments. In this paper, we provide 
a quantitative assessment of the uncertainty associated with the theoretical description of the nuclear 
cross sections, estimating it from the discrepancies between the predictions of different models. 

PACS numbers: 13.15.-|-g, 25.30.Pt 


I. INTRODUCTION 


Large water-Cherenkov detectors have proven to be a 
powerful tool to address a number of outstanding physics 
issues. Ongoing and future research programs at Super- 
and Hyper-Kamiokande [I| include searches of diffuse su¬ 
pernova neutrinos (DSN) Si and sterile neutrinos in 
the beam of the T2K experiment Q. 

These studies require that the backgrounds i related 
to atmospheric-neutrino interactions be estimated with 
a challenging level of accuracy, because the signals are 
expected to be elusive—a few DSN events per year in the 
fiducial volume of the Super-Kamiokande (SK) detector. 


The DSN flux consists of neutrinos and antineutrinos 
of all the flavors, produced in the past core-collapse su¬ 
pernova explosions. As in the scale of the whole Universe 
such events take place approximately every second @ 
and have isotropic distribution, the DSN signal is be¬ 
lieved to be steady in time and uniform in space [ 3 , @ 1 . 
Its measurement would shed light on average supernova 
features impossible to obtain otherwise, the core-collapse 
event rate within our Galaxy being only a few per cen¬ 
tury HP- Therefore, the DSN signal provides a unique 
window on the bulk picture of the entire supernova popu¬ 
lation, the understanding of which is essential to address 
open questions of paramount importance, such as the 
origin of heavy chemical elements and the production of 
cosmic rays [ 11 |. 


The DSN search of the SK experiment is currentl y p er- 
formed in the energy range 13 < < 90 MeV @,3, in 

which the dominant interaction mechanism, sensitiye to 
i>e’s only, is inyerse P decay of free protons in the water 


molecule [T3 |. 

De + p ^ e'^-\-n. (1) 

Until recently [2, neutrons could not be obseryed and 
the only signature of the DSN eyent o was the signal of 
the positron, which in water-Cherenkoy detectors cannot 
be distinguished from the signal of an electron or a 7 ray. 
As a consequence, all processes yielding e+’s, e“’s, or 7 ’s 
in the relevant energy range, generated backgrounds to 
the DSN search. 

At present Q, the upgraded electronics of the SK de¬ 
tector allows neutrons to be detected by the measurement 
of 2.2 MeV 7 rays, resulting from their capture on free 
protons, with efficiency ^18%. The requirement that the 
positron-like signal be followed by the neutron-capture 
signature has led to a sizable background reduction, en¬ 
abling the DSN search to be performed down to energy 
~13 MeV 3 , compared to '^17 MeV of the previous 
analysis 0 . 

Further progress will be possible if the gadolinium dop¬ 
ing program, which may boost the neutron capture effi¬ 
ciency to ~90%, will be successfully completed in the SK 
experiment [l3[4l5l| . Within this scenario, the neutron 
produced in reaction o will be signaled by the 8 MeV 
7 -ray cascade from the capture in the gadolinium nucleus, 
easier to observe than the 2.2 MeV emission from the ab¬ 
sorption on free protons. Efficient neutron detection will 
dramatically reduce many backgrounds 0 , allowing the 
DSN search to be extended further into the low-energy 
region, of crucial importance. However, this is not the 
case for the contribution of mechanisms producing 7 rays 
in coincidence with the neutrons, for which precise esti¬ 
mates of the cross sections and their uncertainties are 
required. 
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In this context, the most significant role is played 
by neutron knockout from oxygen induced by neutral- 
current (NC) quasielastic (QE) scattering of atmospheric 
neutrinos and antineutrinos |16l| . This reaction, occur¬ 
ring at a rate of ~2 eventsper day in the fiducial vol¬ 
ume of Super-Kamiokande [l^ , may yield 7 rays through 
two mechanisms: (i) deexcitation of the residual nucleus 
and (ii) interactions of the knocked-out nucleon with 
the surrounding medium. The associated cross sections 
and their uncertainties have been estimated in Refs. [13 
and [ 1 ^ within the approach based on the impulse ap¬ 
proximation scheme and realistic nuclear spectral func¬ 
tions 0123 . In the following, we will refer to it as the 
SF approach. It is worth noting that the predictions of 
Refs. [13, [l3 have been found consistent with the exper¬ 
imental result of the T2K Collaboration [2l| . 

To gauge the uncertainties associated with the the¬ 
oretical description of the nucleon-knockout cross sec¬ 
tions, in this pap er we compare the results of the SF 
formalism to those obtained from different ap¬ 

proaches, namely relativistic plane-wave impulse approx¬ 
imation (RPWIA), relativistic mean-field (RMF) 
relativistic Green’s functions (RGF) [13, and super¬ 
scaling (SuSA) [ 13,113 ■ For the sake of completeness, we 
also show the results of the relativistic Fermi gas (RFC) 
model of Ref. [13 ■ 

Motivated by the requirements of both the ongoing 
search of the DSN signal in SK and the determination 
of the total active-neutrino flux in the T2K experiment, 
we consider both neutron and proton knockout from oxy¬ 
gen induced by (anti)neutrino NC QE interaction. Note 
that NC interactions are not sensitive to neutrino flavor. 
Hence, oscillations between active neutrinos do not affect 
the rate of NC events detected with a near and a far de¬ 
tector. Should a deficit of NC events in the far detector 
be observed, this would point to oscillations into sterile 
neutrinos. To determine the rate of NC QE interactions 
in the kinematical setup of T2K, it is necessary to mea¬ 
sure the associated 7 -ray production, owing to the ^50% 
contribution of neutron knockout and high Cherenkov 
threshold of protons in water [ 3 ]. 

The different theoretical approaches employed to de¬ 
scribe neutrino- and antineutrino-nucleus interactions are 
outlined in Sec. m while Sec.Hniis devoted to the discus¬ 
sion of the calculated NC QE cross sections for oxygen. 
Finally, in Sec. EYl we summarize our findings and state 
the conclusions. 


II. THEORETICAL APPROACHES 

In scattering processes off nuclei, the internal structure 
of the target is probed with space resolution ^ 27r/|q|, 
where |q| denotes the momentum transfer. Therefore, 
for |q| larger than ~ 27r/fi, d being the average inter¬ 
particle distance in the target, nuclear scattering is ex¬ 
pected to reduce to the incoherent sum of elementary 
processes, involving bound moving nucleons [^ . This is 


the premise underlying the impulse approximation (lA), 
providing the conceptual framework of all approaches dis¬ 
cussed in this work. 


A. Superscaling 

The SuperScaling Approch (SuSA) to neutrino scatter¬ 
ing, proposed in Refs. [^, [^ for CC and NC reactions, 
respectively, exploits the measured electron-nucleus cross 
sections to predict neutrino-nucleus cross sections. 

The analysis of electron scattering data in the QE sec¬ 
tor [ 23 , 113 ] shows that, at large momentum transfer, the 
reduced inclusive cross section (i.e. the nuclear cross 
section divided by the sum of the elementary electron- 
nucleon cross sections) is largely independent of both jqj 
(scaling of the first kind) and the nuclear target (scaling 
of the second kind), when represented as a function of 
the scaling variable ip- 

These properties are clearly observed in the longitu¬ 
dinal channel, whereas violations associated with reac¬ 
tion mechanisms not taken into account within the lA 
scheme—such as inelastic scattering and processes in¬ 
volving meson-exchange currents—show up in the trans¬ 
verse channel, mainly at energy transfer larger than that 
corresponding to single-nucleon knockout. 

On the basis of the above observations, a phenomeno¬ 
logical superscaling function /(i/)) has been extracted 
from electron-scattering data within a fully relativistic 
framework. The function f{ip) embodies the essential 
nuclear dynamics, including both initial- and final-state 
physics. Its most striking features are a pronounced 
asymmetric tail at large energy transfer and a maximum 
~20% lower than the RFC prediction. 

Within SuSA, neutrino-nucleus cross sections are sim¬ 
ply obtained multiplying the function /(V') by the ap¬ 
propriate elementary weak cross sections. Although phe¬ 
nomenological, this approach has several merits: (i) it 
agrees by construction (up to scaling violations) with 
electron scattering data, (ii) it accounts for both kinemat¬ 
ical and dynamical relativistic effects—which are known 
to be significant even at moderate momentum and energy 
transfer—and can therefore be safely used in a broad en¬ 
ergy range, and (iii) owing to scaling of second-kind, it 
allows for a consistent treatment of different target nu¬ 
clei. In addition, it is worth mentioning that the model 
has been successfully extended to higher energies, well 
beyond the QE regime [Slj. 

B. Relativistic Approaches 

The RPWIA, RMF and RGF approaches provide a 
microscopic and fully relativistic description of the scat¬ 
tering process, at both kinematical and dynamical level. 

In all cases, the bound nucleon states are represented 
by four-spinors, obtained from the self-consistent solution 
of the Dirac-Hartree equation derived from a Lagrangian 
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including cr, uj and p mesons within the mean-field ap¬ 
proximation [H, . The state of the outgoing nucleon 

is described by a relativistic (four-spinor) scattering wave 
function. 

The RPWIA approach can be regarded as the simplest 
implementation of the above formalism. Within RPWIA 
all final-state interactions (FSI) between the nucleon in¬ 
teracting with the beam particle and the spectators are 
neglected, and the knocked out particle is described by 
a plane wave. This approximation, while leading to a 
dramatic simplification in the description of the process, 
results in the appearance of significant discrepancies be¬ 
tween the calculated cross sections and scaling functions 
and the data. For example, the RPWIA scaling func¬ 
tion does not exhibit the asymmetry clearly visible in 
the data, thus suggesting that a more realistic descrip¬ 
tion, taking into account the effects of FSI, is needed. 

In this work, we have used two different schemes, re¬ 
ferred to as RMF and RGF. The former makes use of 
relativistic distorted waves obtained with the same rela¬ 
tivistic scalar and vector energy-independent potentials 
used to determine the initial bound states. The results 
of this approach are in excellent agreement with the phe¬ 
nomenological scaling functions extracted from electron- 
nucleus scattering data Moreover, the transverse 

scaling function exhibits an enhancement of ~20%, with 
respect to the longitudinal one [s^ . It has been suggested 
that a similar enhancement may have a si gnifi cant impact 
on the measured neutrino cross sections 1361-1^. 


C. Spectral function 

The SF approach is based on the factorisation ansatz, 
which amounts to writing the nuclear final state as a 
product of a plane wave, describing the motion of the 
knocked out nucleon, and a (A—l)-nucleon state, describ¬ 
ing the recoiling spectator system. Within this scheme, 
the nuclear cross section reduces to the convolution of the 
elementary scattering cross section with the target spec¬ 
tral function, yielding the energy and momentum distri¬ 
bution of the struck nucleon in the initial state. 

The spectral function of Refs. [l^,[13, employed in this 
work, accounts for the shell structure of the oxy gen nu¬ 
cleus, extracted from experimental (e, e'p) data (49|, as 
well as for the contribution of short-range correlations be¬ 
tween nucleons, taken from theoretical calculations [^ . 
Those two components have been consistently combined 
by the authors of Refs, [l^, [2^ within the framework 
of the local-density approximation (LDA). The LDA ap¬ 
proach, based on the tenet that nucleon-nucleon correla¬ 
tions are largely unaffected by surface and shell effects, 
allows one to obtain the correlation contribution for a fi¬ 
nite nucleus from the corresponding results computed for 
uniform nuclear matter at different densities. 

The LDA spectral functions have been successfully 
used in the analysis of inclusive electron-scattering data 
for carbon and oxygen targets at beam energies up to 
few GeV [l^, [l^. Moreover, the LDA momentum dis¬ 
tribution of carbon turns out to be consistent with that 
extracted from (e, e'p) data at large missing energy and 
missing momentum [SiL l52j| . 


In the RGF approach, FSI are included using the rela¬ 
tivistic Green’s function formalism and a complex optical 
potential (dol - l^ . The use of a complex optical potential 
allows for a fully consistent description of the inclusive re¬ 
sponse, as the loss of flux associated with the occurrence 
of inelastic processes is taken into account. Moreover, 
because of the analyticity properties of the optical po¬ 
tent!^ the Goulomb sum rule is fulfilled by construction 

[laiii. 


Like the RMF approach, the RGF provides a good 
description of (e, e') data in the QE sector, reproduc¬ 
ing both the asymmetry of the scaling function and the 
transverse enhancement 4^. The RGF approach is ap¬ 
propriate for an inclusive process where only the final 
lepton is detected. On the other hand, in the NG QE 
scattering, only the final nucleon can be detected and 
the cross section is semi-inclusive in the hadronic sector. 
The RGF may therefore include channels which are not 
present in the experimental NG QE cross sections, but it 
can also recover important contributions which are not 
taken into account by other models based on the lA. For 
instance, in comparison with the MiniBooNE NG QE v 
and V scattering data (dsl . , the RMF generally under¬ 
predicts the experimental cross section, while the RGF 
results are in reasonable agreement with the data (^1^ . 


III. RESULTS 

All calculations presented in this article have been per¬ 
formed u sing the same values of physical constants, given 
in Refs. [Tj, We employ the electromagnetic form 
factors parametrized according to Refs. (s^. [54|. and the 
dipole parametrization of the axial form factor 

^ (i + q7m2)2’ 

with the coupling constant pA = —1.2701 [s^. 

For the axial mass Ma, we apply the value 1.03 GeV 
, in g ood agreement with that obtained in recent anal¬ 
yses [13 of deuteron measurements This choice 

appears to be best suited for the purpose of the present 
analysis, since it allows for a clear separation between 
the dynamics of the elementary interaction vertex and 
nuclear dynamics, determining the reaction mechanism. 

Note that in the studies of Ref. mill, aimed at pro¬ 
viding an estimate of measurable cross sections includ¬ 
ing the contribution of multinucleon processes, the axial 
mass was instead set to the value Ma = 1.2 GeV, de¬ 
termined by the K2K Gollaboration for the oxygen 
target. 
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FIG. 1. (Color online) Cross sections for neutron [(a) and (c)] 
and proton [(b) and (d)] knockout from oxygen induced by 
NC QE interaction of neutrino [(a) and (b)] and antineutrino 
[(c) and (d)[. The RPWIA and SF results are shown by the 
dashed lines and the filled circles, respectively. The RFG 
calculations, represented by dotted lines, are also included, 
for reference. 


FIG. 2. (Color online) Same as in Fig. [T] but for the SuSA 
model (dotted lines), the RMF approach (solid lines), and 
the RGF calculations with the EDAI (long-dashed lines) and 
“democratic” (short-dashed lines) potentials. For compari¬ 
son, we also include the SF results displayed in Fig. [T1 repre¬ 
sented by filled circles. 


The RGF calculations have been carried out using two 
parametrizations of the relativistic optical potential of 
the energy-dependent and A-independent (EDAI) 
model of Ref. and the more recent “democratic” 
phenomenological optical potential of Ref. . The for¬ 
mer is a single-nucleus parametrization, constructed to 
reproduce elastic proton-oxygen data, while the latter 
has been obtained from a global fit to more than two 


hundred data sets of elastic proton-nucleus data for a 
broad range of targets, from helium to lead. 

The ^gO(j^, i^iV) and ^gO(P, PiV) cross sections com¬ 
puted using the theoretical approaches discussed in Sec. 
HI] are displayed in Figs. [T] and [21 

In Fig. (H we show the cross sections obtained from 
the RPWIA, RFG, and SF models. Because these ap¬ 
proaches do not include FSI, comparing their results pro¬ 
vides information on the uncertainty associated with the 
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description of the initial state. However, it should be 
kept in mind that—to the extent to which the sum over 
the undetected hadronic final states spans a complete 
set—the total inclusive cross section in the QE sector is 
unaffected by the inclusion of ESI. In the SE approach 
this property is fulfilled by construction. However, in 
other approaches the inclusion of ESI is found to have 
non-negligible effects (see, e.g, [sM 1^1. 

Eor neutrino scattering, the discrepancy between the 
predictions of the different approaches is 21% (23%) at 
0.6 (1.5) GeV, and turns out to be similar for neutron 
and proton knockout. Eor antineutrino scattering, the 
spread is more significant, reaching 40% (29%) at 0.6 
(1.5) GeV for neutron knockout. This feature is not sur¬ 
prising, since in this case destructive interference between 
the response functions makes the cross section more sen¬ 
sitive to differences in modeling. 

In Eig. [21 we compare the results of the phenomeno¬ 
logical SuSA approach to those obtained from the RME 
and RGE models, in which the effects of ESI are explic¬ 
itly taken into account. In the RME model, real scalar 
and vector relativistic potentials are used, whereas the 
EDAI and “democratic” complex optical potentials are 
employed in the REG approach. Note that, unlike the 
previous ones, the implementation of the SuSA model 
employed to obtain the results of Eig. [2] includes the ef¬ 
fects of Pauli blocking, following the procedure described 
in Ref. (^ . For comparison, we also show the results of 
the SF approach, which are identical to those displayed 
in Fig. HI 

It is apparent that, compared to the RME model, the 
SuSA approach yields sizeably lower cross sections. This 
is likely to be ascribed to the fact that the RME predicts 
an enhancement of the transverse response, as an effect of 
off-shell spinor distortion [^. The transverse enhance¬ 
ment, clearly observed in electron scattering data, is not 
reproduced by the present version of the SuSA model, in 
which the same scaling function is used in the longitudi¬ 
nal and transverse channels. Work aimed at improving 
the SuSA approach to implement this feature is discussed 
in Ref. [6^ . It is interesting to observe that the RME and 
SF approaches, while being based on very different mod¬ 
els of nuclear dynamics, yield remarkably similar results. 
The two curves corresponding to the cross sections ob¬ 
tained from the RGE approach lay significantly above the 
ones corresponding to the RME, SuSA and SF models. In 
addition, they show a sizable sensitivity to the the optical 
potentials, the discrepancy between the results obtained 
using the EDAI and “democratic” parametrizations rang¬ 
ing between ^5-10% at Ej, ^ 2 GeV and ~10”25% in the 
energy region ^0.3-0.5 GeV. 

As pointed out in Ref. [l^, the larger cross sections 
in the RGE may be associated with the redistribution of 
the strength arising from reaction mechanisms other than 
single nucleon knockout, such as rescattering of the out¬ 
going nucleon—possibly leading to the excitation of non- 
nucleonic degrees of freedom—or scattering off a nucleon 
belonging to a correlated pair. These channels, although 


not explicitly included in the RGE, may be phenomeno¬ 
logically taken into account by the imaginary part of the 
optical potential. 


IV. CONCLUSIONS 

We have calculated the ^|0(z^, i/V) and ^gO(P, i'N) 
cross sections in the QE channel, using a variety of theo¬ 
retical approaches extensively validated through compar¬ 
isons to electron scattering data. 

We emphasize that our analysis is not meant to assess 
the validity of the assumptions underlying the different 
approaches. The aim of this article is to quantify the 
uncertainty associated with the theoretical description 
of the nuclear cross sections, comparing the predictions 
of the models outlined in Sec. m 

The spread of the theoretical results depends apprecia¬ 
bly on energy, and turns out to be larger for antineutrino 
cross sections. While for the ^80(i?, Pp) knockout and 
Ej, > 1.23 GeV, the highest cross section is predicted 
by the RGE model with the EDAI potential, for all other 
processes and kinematics, the highest results are obtained 
using the RGE model with the “democratic” potential. 
In all the cases, the SuSA calculations yield the lowest 
cross sections. For example, at E^, = 600 MeV, the SuSA 
(z/, ly'N) and (b, P'N) results are lower than the RGE ones 
by ^37% and ^^47%, respectively. When the energy in¬ 
creases, the differences are somewhat reduced, and at 
Ej, = 1.5 GeV, the SuSA cross sections are ^30-33% 
lower than the RGE ones. 

It is apparent that the broad spread of the results is to 
be mainly ascribed to the large cross sections obtained 
from the RGE approach. The discrepancies between the 
predictions of the SF, RME, and SuSA approaches turn 
out to be much less pronounced, not exceeding ^15% 
over the energy range 0.3 < < 2.0 GeV. 

Note that the calculations discussed in this article have 
been performed using the dipole parametrization of the 
axial form factor. While another dependence [6^-[7lj 
may introduce noneligible effects on the cross sections, an 
estimate of the corresponding uncertainties would require 
new experimental constraints, preferably from neutrino- 
deuteron and antineutrino-proton scattering. 

We recall that uncertainty of the strange axial cou¬ 
pling constant. As = —0.08 ±0.05, translates into uncer¬ 
tainties of the ^gO(z/, vN) and ^80(z/, vN) cross sections, 
shown in Ref. [1^ not to exceed ^6% for neutrinos and 
~8% for antineutrinos. On the other hand, its effect 
on the ^ 80(^1 S'lid ^|0(z/, v) cross sections is less pro¬ 
nounced, because the proton and neutron contributions 
largely cancel [H, 113, ■ 

As a final cautionary remark, we point out that multi¬ 
nucleon emission processes, not taken explicitly into ac¬ 
count in any of the approaches outlined in Sec. im are 
known to provide a non-negligible contribution to the nu¬ 
clear cross sections in the QE channel. Their inclusion 
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